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The creation and detection of atomic polarization is examined theoretically, through the study of 
basic optical-pumping mechanisms and absorption and fluorescence measurements, and the depen- 
dence of these processes on the size of ground- and excited-state hyperfine splittings is determined. 
The consequences of this dependence are studied in more detail for the case of nonlinear magneto- 
optical rotation in the Faraday geometry (an effect requiring the creation and detection of rank-two 
polarization in the ground state) with alkali atoms. Analytic formulas for the optical rotation signal 
under various experimental conditions are presented. 
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I. INTRODUCTION 

Since the pioneering work of Alfred Kastler and Jean 
Brossel in the 1950s jT], atomic polarization created by 
the interaction of light with atoms have been an excit- 
ing topic of research, providing new methods for laser 
spectroscopy and delivering new technologies for practi- 
cal applications, such as narrow-band optical filters |2J. 

Atomic polarization, created in a medium by polarized 
light, can modify the optical response of the medium, af- 
fecting the light field. For example, the absorption of 
light of a particular polarization by atoms in a polar- 
ized state can be reduced (electromagnetically induced 
transparency [3]) or increased (electromagnetically in- 
duced absorption [3]) compared to that for an unpolar- 
ized state. Coherent population trapping [5. is a closely 
related phenomenon, the study of which led to the dis- 
covery of an interesting effect that is also a powerful tool 
for the manipulation of atomic states: coherent popula- 
tion transfer between atomic states, known as STIRAP 
(stimulated Raman adiabatic passage) |6j. "Lasing with- 
out inversion" [7J |S] is another related effect . 

Additional effects are encountered when atoms inter- 
act with coherent light in the presence of a magnetic 
field [5J [TO]. (Reference [TT] discusses a relationship 
between these effects and electromagnetically induced 
absorption.) These magneto-optical effects — especially 
those involving magnetic-field- induced evolution of long- 
lived ground-state polarization — can be used to perform 
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sensitive magnetometery [15]. (These effects are also of- 
ten referred to as "coherence effects" , although this is 
something of a misnomer, as in some cases the effects can 
be described using a basis in which there are no ground- 
state coherences [T3].) 

The atomic polarization responsible for specific effects, 
such as nonlinear magneto-optical rotation (NMOR) , can 
be described in terms of the polarization moments (PM) 
in the multipolc expansion of the density matrix |14IU5j . 
The lowest-rank multipolc moments correspond to popu- 
lation, described by a rank k = tensor, orientation, de- 
scribed by a rank k = 1 tensor, and alignment, described 
by a rank k — 2 tensor. It is these three lowest-rank 
multipolc moments that can directly affect light absorp- 
tion and laser-induced fluorescence [T5J [TB] , and thus can 
be created and detected through single-photon interac- 
tions. An atomic state with total angular momentum F 
can support multipolc moments with rank up to k = 2F 
|14l 115] ; multi-photon interactions and multipole transi- 
tions higher than dipole allow the higher-order moments 
to be created and detected. Magneto-optical techniques 
can be used to selectively address individual high-rank 
multipolcs [17, 18 . Recently, the possibility of using the 
k = 4 hexadecapole moment to improve the characteris- 
tics of atomic magnetometers was studied (see Ref. [T9] 
and references therein). Effects due to the k = 6 hexa- 
contatetrapole moment have also been observed [2TJ1 121) . 

Magneto-optical coherence effects that involve linearly 
polarized light generally require the production and de- 
tection of polarization corresponding to atomic align- 
ment. (There are multi-field, high-light-power effects in 
which alignment is converted to orientation, which is then 
detected [22j |23j H3]; these effects still depend on the 
creation of alignment by the light.) Thus, for ground- 
state coherence effects, the ground state in question must 
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have angular momentum of at least F = 1 in order 
to support a rank- two polarization moment. The al- 
kali atoms K, Rb, and Cs — commonly used for magneto- 
optical experiments — each have ground-state hyperfinc 
sublevels with F > 1. If light is tuned to a suitable tran- 
sition between a ground-state and an excited-state hy- 
perfine sublevel, alignment can be created and detected 
in the ground state. 

The situation changes, however, if the hyperfine struc- 
ture is not resolved. If the hyperfine transitions are 
completely unresolved (as was the case in early work 
that used broadband light sources such as electrode- 
less discharge lamps to excite atoms), then it is the 
fine-structure transition that is effectively excited — the 
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■P1/2) or the D2 line (n 2 S 1/2 



n P 3 / 2 ). In this case, the effects related to the excita- 
tion of a particular hyperfine transition are averaged out 
when all transitions are summed over. Thus the effect of 
the nuclear spin is removed, and the states have effective 
total angular momentum J = 1/2 for the ground state 
and J = 1/2 or 3/2 for the excited state. In this case the 
highest rank multipolc moment that can be supported by 
the ground state is orientation (k = 2 J = 1), and effects 
depending on atomic ground-state alignment will not be 
apparent. 

In practical experiments with alkali atoms in vapor 
cells, even when narrow-band laser excitation is used, the 
hyperfine structure is in general only partially resolved, 
due to Doppler broadening. At room temperature, the 
Doppler widths of the atomic transitions in K, Rb, and 
Cs range from 463 MHz for K to 226 MHz for Cs. The 
ground-state hyperfine splittings, ranging from 462 MHz 
for K to 9.192 GHz for Cs, are on the order of or greater 
than the Doppler widths, while the excited-state hyper- 
fine splittings, ranging from 8 MHz to 1.167 GHz, are 
generally on the order of or smaller than the Doppler 
width. Thus the question arises: how do coherence effects 
depend on the ground- and excited-state hyperfine split- 
ting when the hyperfine structure is neither completely 
resolved nor completely unresolved? 

In Sec. |H] we discuss transitions for which one or the 
other of the excited- or ground-state hfs is completely un- 
resolved. We determine which polarization moments can 
be created in the ground state via single-photon interac- 
tions, and which moments can be detected through their 
influence on light absorption. We find that the two con- 
tributions to the ground-state polarization — absorption 
and polarization transfer through spontaneous decay — 
depend differently on the ground- and excited-state hy- 
perfine structure. 



In Sec. Ill we choose a particular system (the Dl and 



D2 lines of alkali atoms) and investigate the detailed de- 
pendence of NMOR signals on the excited- and ground- 
state hyperfine splitting. We consider three cases: sys- 
tems in which the atomic Doppler distribution can be 
neglected, and systems in which the Doppler distribu- 
tion is broad compared to the natural line width and 
in which the rate of velocity-changing collisions is either 



much slower than or much faster than the ground-state 
polarization relaxation rate. Appendix [A| contains some 
general results used in Sec. |III|and some more details of 
the calculation. 

Throughout the discussion we use the low-light- 
intensity approximation in order to simplify the calcula- 
tions and obtain analytic results. It can be shown, using 
higher-order perturbation theory and numerical calcula- 
tions, that the essential results presented here hold for 
arbitrary light intensity, as well. Previous work that dis- 
cusses the dependence of optical pumping on whether or 
not hfs is resolved includes Refs. [231 [2^j |2"T[ |2"51 12"9] . 



II. TOTALLY UNRESOLVED GROUND- OR 
EXCITED-STATE HYPERFINE STRUCTURE 

In this section, we discuss the creation and detec- 
tion of atomic polarization in systems for which either 
the ground- or excited-state hyperfine structure is un- 
resolved. This section deals with systems that can be 
described using the complete-mixing approximation, i.e., 
the assumption that atomic velocities are completely 
rethermalized in between optical pumping and prob- 
ing. This is the case for experiments using buffer-gas 
or antirelaxation-coated vapor cells, in which atoms un- 
dergo frequent velocity-changing collisions during the 
ground-state polarization lifetime. The consequences of 
the complete mixing approximation are similar to those 
of the broad-line approximation [301 I3U [Ml 123] > which 
takes the spectrum of the pump light to be broader than 
the Doppler width of the ensemble. In the complete- 
mixing case, narrow-band light produces polarization in a 
single velocity group in the Doppler distribution, and the 
polarization is averaged over all velocity groups through 
rcthermalization, while in the broad-line case, the entire 
Doppler distribution is pumped directly. 



A. Depopulation pumping 

We consider an ensemble of atoms with nuclear spin 
J, a ground state with electronic angular momentum J g , 
and an excited state with angular momentum J e . The 
various ground- and excited-state hyperfine levels are la- 
beled by F g and F e , respectively. The atoms are subject 
to weak monochromatic light with complex polarization 
vector e and frequency oj, near-resonant with the atomic 
transition frequency ujj g j e . We assume that the atoms 
undergo collisions that mix different components of the 
Doppler distribution (or, equivalently, use the broad-line 
approximation) . We also neglect coherences between dif- 
ferent ground-state or different excited-state hyperfine 
levels (these coherences will not develop for low light 
power as long as the hyperfine splittings are larger than 
the natural width of the excited state). We first consider 
polarization produced in the ground state due to atoms 
absorbing light and being transferred to the excited state 



3 



< 



~i 1 r 

M 
/ \ 



F g = 3 -» F e 
i i i_ 



F„=2^F e 



-2 2 

Detuning (GHz) 

FIG. 1: Doppler-free (solid line) and Doppler-broadened 
(dashed line) absorption spectra for the S5 Rb D2 line. A 
Maxwellian velocity distribution at room temperature is as- 
sumed. If the incident light frequency is tuned near the center 
of the F g = 2 — *• F e transition group, the condition discussed 
in the text is fulfilled. Namely, the light detuning from each 
resonance frequency is either much less than or much greater 
than the Doppler width. The condition holds somewhat less 
rigorously for light tuned to the center of the F g = 3 — > F e 
transition group. 



(depopulation pumping). The general form of the con- 
tribution to the ground-state density matrix due to this 
effect is given by [26] (see also Refs. [301 El 1321 [S3] for 
the derivation in the context of the broad-line approxi- 
mation) 



d mr e ■ d rn G(u> 



UJ, 



0. (i) 



where m and n are degenerate ground states, r is an 
excited state, u rn is the transition frequency between r 
and n, and G is a function describing the spectral line 
shape. If the natural width of the excited state is much 



smaller than the Doppler width Tjj , G is approximately a 
Gaussian of the Doppler width. For the system described 
above, this takes the form 
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{F e m"\e-d\F g m') G(uj-uj FaFg ). 



(2) 



Now suppose that the light frequency is tuned so that 
it is close, compared with the Doppler width, to an unre- 
solved group of transition frequencies, and far from every 
other transition frequency (Fig. [T]) . We employ the sim- 
plest approximation that G(u) — cof c F g ) takes the same 
value for each transition in the unresolved group, and is 
zero for all other transitions. With these approximations, 
Eq. ([2| becomes 
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where the sum now runs over only those excited states F e 
that connect via one of the unresolved resonant transi- 
tions to the ground state F g in question. (This sum also 
arises in the broad- line approximation.) 

We now investigate which coherences can be created 
in the ground state by the light. As we will see, this 
will determine which polarization moments can be cre- 
ated. Suppose first that the excited-state hfs is en- 
tirely unresolved. Then the sum over \F e m"){F e m"\ 
in Eq. (§ runs over all excited states, so that it is 
equivalent to the identity. We replace this sum with 
the sum over the eigenstates in the uncoupled basis 
J2m"m" \Fm'{ J e m'j) (Irrij J e m'j\. Further, we insert ad- 
ditional sums to expand the ground-state coupled-basis 
eigenstates in terms of the uncoupled basis. We also ex- 
pand e and d in terms of their spherical components. 
Equation ^ becomes 



x (Im I J e mj\d- q >>\Im'jJgm'j)(Irn' I Jgrn'j\F g rri') (4) 
= £(-l) 9 +q (e*) q >e q »(F g m\Im I J g mj)(J g mj\d- q >\ J e mj){J e mj\d- q »\ J g m'j){Im I J g m'j\F g m'), 
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where the inner products (•••]•••) are given by the 
Clebsch-Gordan coefficients, with ( Jy,m s \ Jimi Ji m 2i — 
( J\m\ J2TO2I J3TO3). In the second line we have used the 
fact that the electric-dipole operator is diagonal in the 
nuclear-spin states. 

We now use the Clebsch-Gordan condition m\ + ni2 = 



m3, as well as the related electric-dipole selection rule 

{Jimi\d q \J2m 2 ) = unless mi = m 2 + q, (5) 

to determine which coherences P^mFgin' can ^ e nonzero 
in Eq. Q . Traversing the factors in the last line of Eq. Q 
from left to right, we find that a term in the sum is zero 
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unless 

m = mi + mj, m"j = m.j + q', 
m'j = m"j + q , m = mi + m'j. 

From this we find that 

\m' - m\ = \q' + q"\ < 2. (7) 

We can translate a limit on |Am| directly into a limit 
on the rank k of polarization moments that can be cre- 
ated as follows. The polarization moments are the coeffi- 
cients of the expansion of the density matrix into a sum 
of irreducible tensor operators (a set of operators with 
the rotational symmetries of the spherical harmonics) . A 
polarization moment of rank k has 2k + 1 components 
with projections q — —K, . . . , k, which are related to the 
Zeeman-basis density-matrix elements by 

oc J2 (FmKq\Fm') Pmm ,. (8) 

m,m'=—F 

From Eq. ^ , a ground-state PM p q K ^ with a given value 
of | q | can exist if and only if there is a |Am| = |g| co- 
herence in the ground-state density matrix. A limit on 
\q\ is not by itself a limit on k, because any polarization 
moment with rank k > \q\ can have a component with 
projection q. However, if such a high-rank moment exists, 
we can always find a rotated basis such that the compo- 
nent with projection q in the original basis manifests itself 
as a component with projection k in the rotated basis. 
Because Eq. Q holds for arbitrary light polarization, it 
holds in the rotated basis, so we can conclude that no 
polarization moment pq with rank k greater than the 
limit on |Am| can be created, regardless of the value of 

For the case under consideration, this analysis reveals 
that only polarization moments with k < 2 are present. 
This is a consequence of the fact that we are consid- 
ering the lowest-order contribution to optical pumping 
(namely, second order in the incident light field), so that 
multi-photon effects are not taken into account. A single 
photon is a spin-one particle, so it can support polariza- 
tion moments up to k 7 = 2. For a polarization moment 
of rank k to be created, the unpolarized (rank 0) den- 
sity matrix must be coupled to a rank-K PM by the rank 
k 7 < 2 photon. The triangle condition for tensor prod- 
ucts implies that n < k 7 + < 2. 

An additional condition on |Am| can be found from 
Eq. ([6]), using the fact that mj and m'j are projections of 
the ground-state electronic angular momentum, so that 
their absolute values are less than or equal to J g . From 
the first and last conditions of Eq. ^ we find 

\m' -m\ = \m'j-mj\<2J g . (9) 

Thus the coherences that can be created within a ground- 
state hyperfine level F g are limited to twice the ground- 
state electronic angular momentum J s , even if F g > J g . 



As a consequence, polarization moments in the ground 
state are limited to rank k < 2J g . We can under- 
stand this restriction by examining Eq. Q. Because the 
excited-state hyperfine shifts have been eliminated from 
the expression and the electric dipole operator does not 
act on the nuclear spin space, all traces of the hyperfine 
interaction in the excited state have been removed from 
Eq. Q. This is indicated by the fact that, in the last line 
of the equation, the nuclear spin does not appear in the 
state vectors describing excited states. Thus the excited 
state only couples to the electronic spin of the ground 
state, so that there is no mechanism for coupling two 
ground-state nuclear spin states. This means that any 
polarization moment present in the ground state must 
be supported by the electronic spin only. 

Considering now the case in which the excited-state 
hfs is resolved and the ground-state hfs is unresolved, 
opposite to the case considered so far, we find no similar 
restriction. It is clear from Eq. ^ that the polarization 
produced in a ground-state hyperfine level is independent 
of all of the other ground-state levels — only one ground- 
state level F g appears in the equation. If the excited-state 
hfs is resolved, then likewise only one excited-state level 
F e appears. Thus pumping on a transition F g — > F e pro- 
duces the same polarization in the level F g as pumping 
on a completely isolated F g — > F e transition, regardless 
of any nearby (unresolved) ground-state hyperfine levels. 
Any polarization moment up to rank k — 2F g can be pro- 
duced, subject to the restriction k < 2 in the lowest-order 
approximation. 

In fact, these results can be obtained without the need 
for any calculations. It is clear that if all the hyperfine 
splittings are set to zero, the nuclear spin is effectively 
nonintcracting, and can be ignored. In this case, only 
polarization moments that can be supported by the elec- 
tronic spin J g can be produced in the ground state. In 
particular, if we consider polarization of a given (degener- 
ate) ground-state hyperfine level, we must have k < 2J g . 
If the ground-state hyperfine splitting is increased, this 
conclusion must remain unchanged, because the light 
only couples the ground states to the excited states; to 
lowest order it does not make any difference what is go- 
ing on in the other ground-state hyperfine levels. If the 
excited-state hyperfine splitting is then increased, the 
various F g — > F e hyperfine transitions become isolated; 
for an isolated transition the limit on the ground-state 
polarization moments is k < 2F g . Thus we see that the 
limit k < 2J g on the ground-state polarization moments 
occurs when the excited state hfs is unresolved, and this 
limit does not depend on whether or not the ground-state 
hfs is resolved. 

The total angular momentum F g can be significantly 
larger than J g . For example, Cs has I — 7/2 and 
J g = 1/2, so that the maximum value of F g is 4. Thus po- 
larization moments up to rank eight can be produced in 
the ground state by depopulation pumping if the excited- 
state hfs is resolved, but only up to rank one if it is 
unresolved [28 . To second order in the light field the 
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FIG. 2: (Color online) Excitation with ^-polarized light on the 



(a) F g 
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and (b) F g 



— > F e = 1 transitions 
of a totally resolved J g = 1/2 — + J e = 1/2 transition with 
I = 1/2. Alignment is produced in the F g — 1 hyperfine level 
in both cases. Relative atomic populations are indicated by 
the number of dots displayed above each ground-state level. 
Relative transition strengths are indicated by the widths of 
the arrows — here the transition strengths are all the same. 
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FIG. 3: (Color online) As Fig. [5] but with excited-state hfs 
unresolved; light is resonant with the F g — 1 — * F e transition 
group. No alignment is produced in the F g = 1 ground state. 



ground-state polarization that can be created is limited 
to at most rank two in any case. However, the ques- 
tion of whether rank-two polarization can be created is 
an important one: ground-state alignment is crucial for 
nonlinear magneto-optical effects with linearly polarized 
light, as we discuss in Sec. |III| 

This situation is illustrated for linearly polarized light 
resonant with an alkali Dl line (J g = J e = 1/2) in 
Figs. [2] and [3] We choose 1 — 1/2 for simplicity, and 
the quantization axis is taken along the direction of the 
light polarization. In Fig. [2] the hfs is completely re- 
solved. Part (a) of the figure shows light resonant with 



the F„ 



1 



F e = transition. Atoms are pumped 
out of the \F g = 1, m = Q) sublevel, producing alignment 
in the F g = 1 state. (Linearly polarized light in the 
absence of other fields can only produce even-rank mo- 
ments, and an F = 1 state can only support polarization 
moments up to rank two; therefore, the anisotropy shown 
in Fig. [2] must correspond to alignment.) If light is reso- 
nant with the F g = 1 — > F e = 1 transition, as in part (b), 
the \F g = l,m = ±1) sublevels are depleted, producing 
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FIG. 4: (Color online) Excitation with z-polarized light on 
the F g = 2 — *• F e transition group of a Dl transition with 
unresolved excited-state hfs. The nuclear spin is I = 3/2. 
No alignment is produced in the F g — 2 hyperfine level. 
The width of each arrow represents the relative transition 
strength, which can be obtained from terms of the sum in 
Eq. §. 



alignment with sign opposite to that in Fig. [2^i. This can 
be contrasted with the case in which the excited-state 
hyperfine structure is completely unresolved, shown in 



Here, all the Zeeman sublevels of the F„ 



Fig. 

state are pumped out equally- 
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the F g 
±1) 



1 



F, 



-the \F g = 1, m = 0) sub- 
= transition, and the 



\F g = l,m= ±1) sublevels on the F g = 1 — > F e = 1 tran- 
sition. (The relative pumping rates, which can be found 
from terms of the sum in Eq. (J3|, are all the same.) Thus 
no imbalance is created in the F g = 1 sublevel popula- 
tions, and no polarization is created in this state. 

The same principle is illustrated for nuclear spin I = 
3/2 in Fig. [4] The excited-state hfs is unresolved, and 
light is resonant with the F g — 2 — > F e transitions. In 
this case, the m = ±1 ground-state sublevels are pumped 
on two different transitions. The total transition strength 
connecting each F g = 2 sublevel to the excited state is the 
same, and so no polarization is produced in the F g = 2 
state. 

The conclusions of this section must be modified when 
polarization produced in the ground state by spontaneous 
emission from the excited state is taken into account. We 
now consider the effect of this mechanism on the ground- 
state polarization (Sec. II B). 



B. Excited state and repopulation pumping 

Through second order in the incident light field (first 
order in light intensity) , there is one additional contribu- 
tion to the ground-state polarization besides the one con- 
sidered in Sec. II A that due to atoms being pumped to 
the excited state and then returning to the ground state 
via spontaneous emission (repopulation pumping). We 
first consider polarization produced in the excited state. 
The general form of the excited-state density matrix is 



^ e • d rk e* ■ d ks G(uj - u r k), 



(10) 
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where r and s are excited states and A: is a ground state. 
Comparing this expression to the formula for ground- 
state depopulation pumping [Eq. (JTJ] , we find that, as one 
would expect, the roles of the ground-state and excited 
state have been reversed. This means that the results of 



Sec. II A with F g and F e interchanged, can be applied 
to the excited-state polarization. In this case, there is 
a limit k < 2J e on the polarization moments that can 
be produced in the excited state, that occurs only when 
the ground state hfs is unresolved. The restriction does 
not depend on whether or not the excited-state hfs is 
resolved. There is the additional limit k < 2 for low light 
power. 

When the polarized atoms in the exited state decay 
due to spontaneous emission, the polarization can be 
transferred to the ground state. This contribution to the 
ground-state density matrix is given by |26] 



(repop) 



,(>r 



with p rs as given above. The fact that this formula has 
no reference to individual transition frequencies leads us 
to expect that the polarization transfer should be inde- 
pendent of the hypcrfine splittings. Indeed, writing this 
expression out for the case under consideration gives 



(repop) 
yF n m,F„; 



<xY^(-l) p (F g m\d p \F e m")(F e m"\p\F e m'' 



(F e m'"\d. p \F g m'} 



(12) 



and the only restriction to be obtained is to' — m = 
to'" — to" (excited-state Am equals ground-state Am) 
[52l 155] . (Transforming to the uncoupled basis does not 
result in any additional limits.) In other words, if the po- 
larization moment can be supported in the ground state, 
it can be transferred from the excited state via sponta- 
neous emission. 

Combining these results, we see that there is a simi- 
lar restriction on polarization created in the ground state 
by repopulation pumping as the one on polarization cre- 
ated by depopulation pumping. However, the restric- 
tion occurs in the opposite case. When the ground state 
is unresolved the polarization produced by repopulation 
pumping must have n < 2J e . This limit does not depend 
on whether the excited-state hfs is resolved. 

We now illustrate the foregoing for a system with 
J g = J e =1=1/2 pumped with linearly polarized light. 
In Fig. [5] both the ground- and excited-state hfs is re- 
solved, and light is tuned to the F g = 1 — > F e = 1 tran- 
sition. In part (a) of the figure, the pump light produces 
polarization in the F e = 1 excited state. In part (b) the 
excited atoms spontaneously decay. This creates polar- 
ization in the F g = 1 ground state, because more atoms 
are transferred to the \F g = 1,to = 0) sublevel than to 
the \F g = l,m = ±1) sublevcls. (In this and the follow- 
ing two figures, we do not show the atoms that decay to 
the F g = state.) Figure [6] is the same but with light 
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FIG. 5: (Color online) Level diagram for an Dl transition 
with resolved hfs for an atom with 1 = 1/2 showing (a) optical 
excitation and (b) spontaneous decay with linearly polarized 



light resonant with the F g 



1 transition. The 



branching ratio for each allowed decay is the same, leading to 
(11) an excess of atoms in the \F g = 1, m = 0) sublevels over the 
populations of the \F g = l,m = ±1} sublevels by a ratio of 
2:1. 
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(Color online) As Fig. [5] but with light tuned to the 
F g = — > F e = 1 transition. In this case an excess of atoms 
results in the m = ±1 states, so that the polarization has the 
opposite sign as that in Fig. [5] 



tuned to the F g = — > F e = 1 transition; polarization is 
also created in the F g = 1 ground state in this case. 

In Fig. [7] the ground-state hfs is now unresolved, while 
the excited-state hfs remains resolved. In this case, both 
ground-state hyperfine levels are pumped by the light, 
and equal populations are produced in the sublevels of 
the F e = 1 state, as shown in part (a) of the figure. As 
seen in part (b), the excited-state atoms spontaneously 
decay in equal numbers to the F g = 1 sublevels, so that 
no polarization is produced in the F g = 1 state. 

Note that in the opposite case, with unresolved 
excited-state hfs and resolved ground-state hfs, sponta- 
neous decay is not prevented from producing polariza- 
tion in the F g = 1 ground state. In this case, atoms are 
pumped into the F e = state, as well as the F e = 1 
state, as shown in Fig. [sJl. Since the \F e = 0, m = 0) 
state decays isotropically, the decay from this state does 
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(a) 

F e = \ 
F, = 



(b) 




F s = 



-1 m = 




FIG. 7: (Color online) As Fig. [5] but with unresolved ground- 
state hfs; light is tuned to the F g — > F e = 1 transition group. 
The contributions to the ground-state polarization illustrated 
in Figs. [5] and [6] cancel, so that no ground-state polarization 
is produced. 



(a) 

F = l . 
F = o 



(b) 




not cancel out the polarization created by decay from the 
F g = 1 state (Fig. [83). Thus we see that it is the ground- 
state hfs, and not the excited-state hfs, that needs to be 
resolved in order for polarization to be produced in the 
ground state due to spontaneous decay. 

To summarize the results obtained so far, to lowest 
order in the excitation light, polarization can be either 
produced in the ground state directly through absorp- 
tion, or transferred to the ground state by spontaneous 
emission. To this order, polarization moments due to 
both of these mechanisms must have rank k < 2. In 
addition, if the excited-state hfs is unresolved, there is 
a limit k < 2J g on the ground-state polarization due to 
depopulation, but no additional limit on the polarization 
due to repopulation. On the other hand, if the ground- 
state hfs is unresolved, there is a limit k < 2 J e on the 
ground-state polarization due to repopulation, but no ad- 
ditional limit on polarization due to depopulation. Thus, 
unless both the excited-state and ground-state hyperfinc 
structure is unresolved, one or the other of the mecha- 
nisms is capable of producing polarization of all ranks 
k < 2. 



C. Light absorption 



F g = 

m = - 1 m = m = 1 m = - 1 m = m = 1 

FIG. 8: As Fig.[5]but with unresolved excited-state hfs; light 
is tuned to the F g = 1 — > F e transition group. Three de- 
cay channels transfer atoms to the \F g = l,m — 0} sublevels, 
while the \F g — l,m — ±1} sublevels are each fed by two de- 
cay channels. Since all the branching ratios are the same, the 
resulting population imbalance is 3:2. 



The absorption A of a weak probe light beam is given 
in terms of the ground-state density matrix by [23] 

■A oc y^e- d rm p mn e* ■ d nr G(uj - Lu rm ), (13) 

mnr 

or 



A(x^2 (Fem\e ■ d\F g m')(F g m'\p\F g m")(F g m"\d ■ e\F e m)G(uj - uj FcFg ), (14) 

where all quantities are as defined above. Using the approximation, as in Sees . |II A] and [TTB| that the light is resonant 
with an unresolved transition group and far detuned from all other transitions, this formula reduces to 

A<*Y^ ( p e m \e ■ d\F g m')(F g m'\p\F g m")(F g m"\d ■ e*\F e m), (15) 



where the sum over F g and F e includes only those com- 
binations that are in the unresolved resonant transition 
group. 

We now investigate the dependence of the absorption 
on the ground-state polarization in various cases. Con- 
sider the case in which the ground-state hfs is completely 
resolved, and the excited-state structure is unresolved. 



The light is tuned to a unresolved transition group con- 
sisting of transitions between one ground-state hyperfinc 
level F„ and all of the excited-state levels. The sum in 



Eq. ( 15 1 over the excited states is then a closure relation, 



and can be replaced with a sum over any complete basis 
for the excited state, in particular, the uncoupled basis. 
We also insert closure relations to expand the ground 
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states (F g m"\ and \F g m') in the uncoupled basis. We obtain 



A<xY^{- 1 ) 9 ' +q ''e q '(e*) q '>{ImiJ e mj\d- q ,\ImU g m'j){Im! I J g ri 

x (Im'j J g m'j\d^ q "\ImjJ e mj) 

= Y.(-v 9 ' +9 '' e *'( e *^''( J e m j\ d -<i'\ J 9 m 'j^ Imi ^ 

(16) 



where only one nuclear-spin summation variable remains 
in the last line. The dipole matrix element selection rules 
and Clebsch-Gordan conditions require that 



m = ni[ + rrtj, m j = mj + q , 
m" = mi + m "ji m J = m "j + ? i 



(17) 



must be satisfied in order for a term in the sum to con- 
tribute to the absorption. These conditions can be com- 
bined to yield \m' - m"\ = \q' + q"\ < 2. Thus only co- 
herences with |Am| < 2 (and polarization moments with 
k < 2) can affect the lowest-order absorption signal. The 
reason for this is analogous to the reason that polariza- 
tion moments of maximum rank two can be created with 
a lowest-order interaction with the light. Absorption oc- 
curs when an atom is transferred to the excited state, 
i.e., when population (rank zero polarization) is created 
in the excited state. Thus, to be observed in the signal, 
a ground-state atomic PM must be coupled to a k = 
excited-state PM by a spin-one photon, which can sup- 
port polarization moments up to rank two. The triangle 
condition for tensor products then implies that the rank 
of the atomic polarization moment must be no greater 



than two. 

Another restriction on the coherences that can affect 



absorption can be found from Eq. ( 17 1 by using the fact 
that I Til j < J„ and \m" T \ < J„. We find 



< 2J„ 



(18) 



In other words, only polarization moments with k < 2 J g 
can affect the absorption signal, regardless of the value 
of F g . Evidently, it is the excited-state hfs that deter- 
mines which ground-state polarization moments can be 
detected in absorption, whether or not the ground-state 
hfs is resolved. 

Considering the case in which both the excited- and 
ground-state hfs is entirely unresolved can lend some in- 
sight into this result. In this case, every combination of 
F g and F e enters in the sum in Eq. (15). If the ground- 



state hyperfine splitting is sent to zero, the sum must be 
extended to include matrix elements of p between differ- 
ent hyperfine levels. This means that all of the sums in 



Eq. ( 15 1 can be replaced with sums over uncoupled basis 



states, giving 



A oc \ (—l) q+q e q i(e*) q ''(ImjJ e mj\d- q i\Im'jJgm'j}(Im l jJgm l j\p\Im'lJ g m'j}(ImjJ g mj\d 

, „ (19) 
= ^(-l) q +q e q >(e*) q ii(J e mj\d^ ql \Jgm'j){Im I Jgm'j\p\Im I J g 7n'j){Jgm'j\d^ q n\J e 7nj). 



Since the hyperfine interaction has been effectively elim- 
inated, the absorption no longer depends on the nuclear 
spin: the complete density matrix does not enter, but 
rather the reduced density matrix 

PrrJjm'j ~ Pm I m' J ,m I m" J (20) 

that is averaged over the nuclear spin mj. The reduced 
density matrix can only support polarization moments up 
to rank k = 2J g , so any PM in p with higher rank cannot 
affect the absorption. Considering a density matrix that 
is nonzero only within one ground-state hyperfine level 



F g , we see that polarization moments with rank greater 
than two will not contribute to the signal. Since the other 
ground-state hyperfine levels are unoccupied, it makes no 
difference what the ground-state hyperfine splitting is, so 
we regain the result that, even if the ground-state hfs is 
resolved, only polarization moments with k < 2J g can 
affect the absorption of light if the excited-state hfs is 
unresolved. 

There is no corresponding restriction on the polar- 
ization moments that can affect absorption when the 
ground-state hfs is unresolved and the excited-state hfs is 
resolved. Indeed, we can consider the case in which only 
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F e = \ 



F e = 



F g = 



(a) (b) 

m = — 1 hi = m = 1 m = - 1 m = m = 1 



FIG. 9: (Color online) Dl transition for an atom with 1 = 1/2 
subject to linearly polarized light resonant with the F g = 1 — * 
F e = 1 transition. In part (a) the F g = 1 ground state is 
unpolarized and there is light absorption. In part (b) the 
F g = 1 state has the same total population, but is aligned, 
and there is no absorption. 



one ground-state hyperfinc level F g is populated: the ab- 
sorption is then exactly as if the transition F g — ► F e were 
completely isolated. For such an isolated transition, the 
only limit on detectable polarization moments is n < 2 
for the low-power case. 

As in the previous subsections, we illustrate this result 
for a Dl transition for an atom with 7=1/2 subject to 
linearly polarized light. In Fig. [9] both the ground- and 
excited-state hfs is resolved, and the light is resonant 
with the F g = 1 — > F e = 1 transition. In part (a) there 
is no polarization in the F g = 1 ground state: atoms are 
equally distributed among the Zeeman sublevels. Light is 



absorbed by atoms in the \F„ = 1, 



±1) sublevels. In 



part (b) there are the same total number of atoms in the 
F g = 1 state, but they are collected in the to — sublevel. 
The population is the same, but the F g — 1 state now 
also has alignment. In this particular case there is no 
absorption, because the atoms are all in the to — dark 
state. Thus, in this situation, the rank-two polarization 
moment has a strong effect on the absorption signal. 



Figure 10 shows the same system, but with unresolved 
excited-state hfs. In this case there is no dark state; all 
of the atoms interact with the light. The distribution of 
the atoms among the Zeeman sublevels does not affect 
the light absorption, and so the rank two polarization 
moment is not detectable in the absorption signal. 



D. Fluorescence 

Finally, we consider which excited-state polarization 
moments can be observed in fluorescence. Assuming 
broadband detection, the intensity of fluorescence into a 
particular polarization e is given in terms of the excited- 
state density matrix by 



X oc ^ 6 • d mr p rs G • d 5m . 



(21) 



(a) 



(b) 



F e = l 



m = - 1 m = m=l m = -l m = m=l 



F, = 



FIG. 10: (Color online) As Fig. [5J but with unresolved 
excited-state hfs. In this case there is no difference in the 
absorption seen for an (a) unpolarized and (b) aligned F g = 1 
ground state. 



Because the sums in r and s go over all excited states, 
and to runs over all ground states, we can write Eq. (21 1 
for our case in terms of the uncoupled-basis states 
gives 



This 



1 oc ^(-1) 9 '+ 9 " (e*) q , e q ,,{Im 1 J g mj\d- g ,\Im' I J e m'j){Im' I J e m' J \p\ImU^^^^^ (22) 

I 



resulting in the restrictions 

m" = -q" + mj, to j = -q 1 + m'j, 
mi = to j = Til j , 

on the terms that can contribute to the fluorescence. This 
indicates that the nuclear polarization cannot affect the 
fluorescence signal, and so only the electronic excited- 
state polarization of rank k < 2 J e can be observed. In 
addition, only coherences with \m'j — m'j\ = \q' + q"\ < 2 



can be observed. This rule has appeared earlier as a 
consequence of the low-light-power assumption; because 
spontaneous decay is not induced by an incident light 
field, in this case the rule is exact. This means that no 
matter the value of J e , and what polarization moments 
exist in the excited state, only polarization of rank k < 2 
can be observed in fluorescence. 
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unresolved limit on k 



Ground-state pol. (depop.) 


excited 


hfs 


< 


2J g 


Ground-state pol. (repop.) 


ground 


hfs 


< 


2Je 


Excited-state pol. 


ground 


hfs 


< 


2J e 


Absorption 


excited 


hfs 


< 


2J g 


Fluorescence 






< 


2J e 



TABLE I: Summary of the results of this section. For each 
quantity, the restriction on the rank k of the polarization 
that can be created or detected is given in the third column. 
The restriction holds under conditions (ground- or excited- 
state hfs unresolved) described in the second column. For 
fluorescence with broad-band detection the restriction holds 
regardless of whether the hfs is resolved. 



E. Summary 

In this section, we have shown that, when the ground- 
or excited-state hfs is unresolved, there are restrictions 
on the rank of the polarization moments that can be 
created or detected by light. Some of these restrictions 
may at first seem counter-intuitive, but they can be ob- 
tained from very basic considerations. For example, the 
two facts that nuclear spin can be ignored if the hfs is 
completely unresolved and that lowest-order depopula- 
tion pumping of a given hypcrfinc level does not depend 
on ground-state hyperfine splitting lead directly to the 
result that polarization moments produced by depopula- 
tion pumping are subject to a limit of k < 2J g when the 
excited-state hfs is unresolved. Various processes of cre- 
ation and detection of polarization are subject to differ- 
ent restrictions (Table]!]). In particular, the two processes 
that can create ground-state polarization — depopulation 
and repopulation pumping — are subject to restrictions 
under different conditions. Consequently, unless the hfs 
is entirely unresolved, there is always a mechanism for 
producing polarization limited in rank only by the total 
angular momentum, rather than the electronic angular 
momentum. 



III. PARTIALLY RESOLVED HYPERFINE 
STRUCTURE: NONLINEAR 
MAGNETO-OPTICAL EFFECTS 

Now let us examine the more general case of partially 
resolved hyperfine transitions. For this study, we will 
look at the quantitative dependence on hyperfine split- 
ting of nonlinear optical rotation — rotation of light po- 
larization due to interaction with a J g — > J e transition 
group in the presence of a magnetic field. In this case, 
the effect of ground-state atomic polarization is brought 
into starker relief: in the experimental situation that we 
consider, both the creation and detection of ground-state 
polarization is required in order to see any signal whatso- 
ever. When linearly polarized light is used, as is supposed 
here, the lowest-order effect depends on rank-two atomic 
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Magnetic 
field 









, - r 






- — — - ^ 






























/ ~ ~ ' ' 

\y- " ~~ 1 

T Transit effect / 



Laser beam 



Antirelaxation-coated vapor cell 

FIG. 11: The transit and wall optical-rotation effects. Both 
effects occur in an anti-relaxation-coated vapor cell, and can 
be distinguished by their greatly different magnetic resonance 
widths. 



alignment. Thus, for the alkali atoms, the question of the 
dependence of the effect on hyperfine structure arises, 
because, as discussed in the previous section, both the 
creation and the detection of alignment in the J g = 1/2 
ground state can be suppressed due to unresolved hfs. (In 
fact, a higher-order effect can occur wherein alignment is 
created, the alignment is converted to orientation, and 
the orientation is detected [22j |23j [24]. However, the 
conversion of alignment to orientation is an effect of ten- 
sor ac Stark shifts, which can be shown by arguments 
similar to those in Sec. [TT] to suffer suppression due to 
unresolved hfs in the same way as does the direct detec- 
tion of alignment.) 

In the Faraday geometry, linearly polarized light prop- 
agates in the direction of an applied magnetic field, and 
the rotation of the light polarization direction is mea- 
sured. A number of magneto-optical effects can con- 
tribute to the optical rotation, including the linear Fara- 
day effect, the Bennett-structure effect, and various ef- 
fects depending on atomic polarization ("coherence ef- 
fects") Here we are concerned with optical rotation 
due to several different forms of the ground-state coher- 
ence effect, in which the atomic velocities are treated 
in three different ways. First we consider the atoms to 
have no velocity spread, and analyze the Doppler-free 
"transit effect", as for an atomic beam with negligible 
transverse velocity distribution |34j . We then consider 
the case in which atoms have a Maxwellian distribution, 
but do not change their velocities in between pumping 
and probing — this corresponds to the transit effect for 
buffer-gas- free, dilute atomic vapors [T3]. Finally, we 
treat the case in which atoms undergo velocity-changing 
collisions between pumping and probing, as for buffer-gas 
cells [35] or the wall-induced Ramsey effect ( "wall effect" ) 
in antirelaxation-coated vapor cells |36j . Figure 11 illus- 

We 



trates the transit and wall effects in a vapor cell. 
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examine the dependence of these effects on the size of 
the hyperfine splittings as they vary from much smaller 
than the natural width to much greater than the Doppler 
width. 

Throughout this section we consider formulas for the 
optical rotation signal valid to lowest order in light power, 
under the assumption that the ground-state relaxation 
rate 7 is much smaller than both the excited-state natu- 
ral width r and the hyperfine splittings. For the Doppler- 
free case a single analytic formula can be applied to both 
resolved and unresolved hfs (i.e., no assumption need 
be made about the relative size of the hyperfine split- 
tings and the natural width) . For the Doppler-broadened 
cases, analytic results can be obtained in various limits, 
which together describe the signal over the entire range 
of hyperfine splittings. 

We first focus on the simplest case: the Dl line 
(J g = J e = 1/2) for an atom with I = 1/2. This is a 
somewhat special case, because one of the two ground- 
state hyperfine levels has F g = 0, and consequently can 
neither support atomic alignment nor produce optical 
rotation. We then consider the differences that arise 
when considering higher nuclear spin and also the D2 
line (J g — 1/2 and J e = 3/2). Finally, results for the 
"real-world" alkali atoms commonly used in experiments 
are presented. Some details of the calculation and gen- 
eral formulas for arbitrary J g , J e , and / are presented in 
Appendix [A] These formulas are generalizations of those 
given in Ref. [T3]; related earlier work includes that of 

Refs. nana. 



A. Doppler-free transit effect 

We consider nonlinear Faraday rotation on a J g — > J e 
atomic transition for an atom with nuclear spin /. We 
can limit our attention to the ground-state coherence ef- 
fects by using a "three-stage" model for Faraday rotation 
1 1 3 j . in which optical pumping, atomic precession, and 
optical probing take place sequentially, and the light and 
magnetic fields are never present at the same time. In 
this case, the linear and Bennett-structure effects, which 
require the simultaneous application of light and mag- 
netic fields, do not occur. Such a model can be realized 
in an atomic beam experiment, but it is also a good ap- 
proximation to a vapor cell experiment that uses low light 
power and small enough magnetic fields so that the co- 
herence effects are dominant. 

The calculation is performed using second order per- 
turbation theory in the basis of the polarizatio n mo ments 



field is turned off, and the light field is applied once more 
to probe the atomic polarization. The nonlinear optical 
rotation is found to lowest order in the probe-light Rabi 
frequency (Appendix A 2). 

Because the magnetic field is neglected during the opti- 
cal pumping stage, the atomic ground-state polarization 
that is produced in this stage is entirely along the light 
polarization direction, i.e., it has polarization component 
q = 0. Since linearly polarized light has a preferred axis, 
but no preferred direction, it can not, in the absence 
of other fields, produce atomic polarization with a pre- 
ferred direction, i.e., polarization with odd rank k. Also, 
we have seen in Sec. [Tljthat, to lowest order in the light 
power, optical pumping cannot produce polarization mo- 
ments with k > 2. Thus the only ground-state polariza- 
tion moment with rank greater than zero that is produced 
at lowest order has n = 2 and q = 0. We first consider 
the Dl line (J g = J e = 1/2) for an atom with I = 1/2. 
In this case, the only ground-state hyperfine level that 
can support the p^ 20 ^(F g F g ) moment has F g = 1. (Due 
to the assumption that the hyperfine splittings are much 
greater than the ground-state relaxation rate, we can ig- 
nore ground-state hyperfine coherences throughout the 



discussion.) From Eq. (AlOl, the value of this moment 
is found to be 



p( 2 °)(ll)= i ^([L(^ 1 )-L« 1 )] 



+ |[L( W i,o)-LKx)]), 



(24) 



where k s — ( J g \\d\\ J £ ) 2 £q /(T-f) is the reduced optical- 
pumping saturation parameter (£0 is the optical electric 
field amplitude), R is the branching ratio for the tran- 
sition J e — > J g , and ^i? e F g i s ^ ne transition frequency 
between excited-state and ground-state hyperfine levels 
in the frame "rotating" at the Doppler-shifted light fre- 
quency lu: d>'p e p = ^FcFg — cij + k • v, where ^F^F g is the 
transition frequency in the lab frame, u> is the light fre- 
quency, k is the wave vector, and v is the atomic velocity. 



We also write 



F, F D 



= — Ai? c i? s + k • v, where Ap e F is 



the light detuning from resonance. We have defined the 
Lorcntzian line profile 



L(J) = 



(r/2) 2 



(r/2)* 



(25) 



Equation (24 1 is written as the sum of two terms, each 



surrounded by square brackets. The first term is the con- 
tribution to the polarization due to depopulation pump- 
ing discussed in Sec. |II A| This term is itself a sum of 
p( K i) (F1F2) of the density matrix (Appendix Al). The contributions due to pumping on the F g — 1 — > F e = 

transition and the F ' g = 1 — > F e = 1 transition. These 



three stages of the calculation are as follows. In stage 
(a), a re-directed light beam linearly polarized along z is 
applied, and we calculate optical pumping through sec- 
ond order in the optical Rabi frequency. In stage (b), 
the light field is removed, and a x-directed magnetic field 
is applied. We calculate the effect of this field on the 
atomic polarization. Finally, in stage (c), the magnetic 



two contributions are of opposite sign, as illustrated in 
Fig. [2] Pumping on cither transition produces alignment 
in the F g — 1 ground state; the sign of the corresponding 
polarization moment depends on whether there is more 
population in the m = sublevel or the m = ±1 sub- 
levels. We saw in the discussion of Sec. Ill Al that when 
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the excited state hfs is unresolved, polarization with rank 
K > 2J g cannot be created by depopulation pumping 
(Fig. [3]). We see here that as w ,i approaches wi t i, i.e., 
as the excited-state hyperfine splitting goes to zero, the 
contributions from the two transitions cancel and this 
term goes to zero. For the Doppler-broadened atomic 
ensemble discussed in Sec. [TTJ the hfs was considered un- 
resolved when the hyperfine splittings were smaller than 
the Doppler width. Since Eq. (24 1 describes a single ve- 



locity group, the relevant width here is the natural width 

r. 



The second term of Eq. (24 1 is the contribution to the 



ground-state polarization due to repopulation pumping 
discussed in Sec. |IIB| This term is also composed of two 
contributions of opposite sign: one due to pumping on 
the F g = 1 — > F e = 1 transition and one due to pump- 
ing on the F g = — > F e — 1 transition. The two con- 
tributions are illustrated in Figs. [5] and [6] wh ich show 
the origin of the opposite signs. In Sec. [II B] we found 
that depopulation pumping cannot create polarization 
moments with rank k > 2J e when the ground-state hfs 
is unresolved (Fig. [7]). We see here that this term of 
Eq. (24) goes to zero when 0*1,0 approaches i.e., as 



the ground-state hyperfine splitting goes to zero. 

In the second and third stages of the model of the co- 
herence effect, the ground-state polarization precesses in 
a magnetic field and is probed by light with the same po- 
larization as the pump light considered in the first stage. 
From Eq. ( A22 1 we find that the normalized optical ro- 



tation dip per path length dt is proportional to the po- 
larization produced in the first stage and is given by 



dip 



jVf[i«i)-iKi)]ariP (20) (ll), (26) 



where 



(7/2)Aj 



(7/2) 5 



1 a 



(27) 




FIG. 12: Illustration of the rotating polarizer model for opti- 
cal rotation, (a) Optical pumping on a F g — 1 — > F e = or 
F g — 1 — > F e — 1 transition causes the medium to act as a 
polarizing filter with transmission axis along the input light 
polarization e. (b) When the transmission axis rotates due 
to Larmor precession, the output light polarization e' follows 
the transmission axis and so rotates in the same sense as the 
Larmor precession, (c) If polarization produced by pump- 
ing on one transition is probed on the other, the polarization 
functions as a polarizing filter with transmission axis perpen- 
dicular to the input light polarization. (Attenuation of the 
light beam is not indicated.) (d) When the medium polar- 
ization rotates, the output light polarization tends to rotate 
towards the transmission axis, in the opposite sense to the 
Larmor precession in this case. 



is the magnetic-resonance line-shape parameter, with 
Qp = gF g fiBB the Larmor frequency for the ground- 
state hyperfine level F g (gp g is the Lande factor for the 
ground state F gi and [Ib is the Bohr magneton), and 



1 dl 
l~d£ 



2tt (2J fl + l) 
RnX 2 (2J e + 1) 



(28) 



is the unsaturated resonant absorption length assuming 
totally unresolved hyperfine structure, where X is the 
light intensity, n is the atomic density, and A is the light 
wavelength. The branching ratio R enters here because 
it factors into the transition strength. 

The contributions to the optical rotation signal from 
the F g = 1 — > F e = transition and the F g = 1 — > F e = 1 
transition have opposite signs. To understand this, it is 
helpful to think of the optically polarized medium as a 
polarizing filter |13j . When pumping on a 1 — > or 
1 — > 1 transition, the medium is pumped into a dark 



(non-absorbing) state for that transition (Fig. ]2j), corre- 
sponding to a polarizing filter with its transmission axis 
along the input light polarization axis e [Fig. 12 'a)]. The 
Larmor precession induced by the magnetic field causes 
the transmission axis of the filter to rotate, so that it is 
no longer along e. This in turn causes the output light 
polarization axis e' to rotate. The polarization of light 
passing through a polarizing filter tends to rotate towards 
the transmission axis, so that in this case the optical 
rotation is in the same sense as the Larmor precession 
[Fig. \L2^b)}. Now, compare the polarization produced 
when pumping on a 1 ^ or 1 -> 1 transition, as shown 
in Fig. [2] We see that the dark state for each transition is 
a bright (absorbing) state for the other. This means that 
if we choose one or the other of these states, it will func- 
tion as just described for one of the transitions, but will 
function as a polarizing filter with its transmission axis 
perpendicular to e for the other transition [Fig. |12[c)]. 
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FIG. 13: Dependence on light detuning from the F g — 1 — » 
_F e = 1 transition of (top) the components of ground-state 
alignment d ue t o depopulation (dashed) and repopulation 
(solid) [Eq. ( 24 1] and (bottom) optical rotation for a given 
amount of alignment [Eq. (26)]. Gray vertical lines show 
F g — > F e transition resonance frequencies. Parameter values 
in units of F are 7 <C 1, A g — 10, A e = 5. 



When the axis of the filter rotates in this case, the fact 
that the output light polarization tends to rotate towards 
the transmission axis, means that here the optical rota- 
tion is in the other direction, in the opposite sense to 
the Larmor precession [Fig. [l2"|d)]. In other words, for 
a particular sign of the rank-two polarization moment, 
the optical rotation will have one sign when probed on 
one transition, and the opposite sign when probed on 
the other, as indicated by Eq. (26 1. Because the observa- 



tion of optical rotation requires the detection of rank-two 
polarization moments, we might expect, analogously to 
the discussion in Sec. 



IIC 



that it is suppressed when 
the excited-state hyperfine splitting goes to zero. Equa- 
tion ( 26 ) shows that the two contributions indeed cancel 
when wq 1 approaches u>i 1. 



Equation ( 26 ) and the two components of Eq. ( 24 ) are 



plotted as a function of light detuning from the F g — 
1 ^ F e ~ 1 transition in Fig. |13| for particular values 
of the ground- and excited-state hyperfine coefficients A g 
and A e . (For J = I = 1/2, the hyperfine coefficient A is 
equal to the splitting between the two hyperfine levels.) 
Here and below numerical values of frequencies are given 
in units of T. As discussed above, each spectrum consists 
of two peaks of equal magnitude and opposite sign. For 
the spectrum of alignment due to depopulation and the 
spectrum of rotation for a given amount of alignment, 
the peaks are separated by the excited-state hyperfine 



FIG. 14: Spectra of normalized optical rotation 
£o/(n a xi)(dip/d£) for the Doppler-free transit effect. Left 
column: components due to polarization produced by de- 
population (dashed) and repopulation (solid); right column: 
total signal. Parameter values in units of F are 7 <C 1, 
A q > l,A e . 



splitting, so that they cancel as this splitting goes to 
zero. For the spectrum of alignment due to repopulation, 
the peaks are separated by the ground-state hyperfine 
splitting; they cancel as the ground-state splitting goes 
to zero. 

In this subsection we are analyzing a Doppler-free sys- 
tem, i.e., we assume that the atoms all have the same ve- 
locity, which we take to be zero for simplicity. Then the 
observed optical rotation signal is found by simply substi- 
tuting Eq. ( 24 1 into Eq. ( 26 ) . We first consider the case 
in which the ground-state hfs is well resolved. The rota- 



tion signal is plotted in Fig. 14 for large ground-state hy- 
perfine splitting and various excited-state splittings A e . 
The components of the rotation signal due to depopula- 
tion (dashed) and repopulation (solid) are plotted in the 
left-hand column, and the total signal is plotted on the 
right. As the previous discussion indicates, the rotation 
signal decreases as the excited-state hyperfine splitting 
A e becomes smaller, with the component due to depop- 
ulation decreasing faster than the component due to re- 
population. This is also seen in Fig. |15| which shows the 
maximum magnitude of the rotation spectrum as a func- 
tion of A e (for each value of A e , the signal is optimized 
with respect to detuning). Thus, for small splittings, the 
component due to repopulation dominates. To lowest 
order in A e , the signal is given by 



dip = A e K,a: 1 J?(r/2) 4 A 1 
dl 144[(r/2) 2 + A2] 3 ' 



(29) 
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FIG. 15: Maximum of the spectrum of the Doppler-free non- 
linear magneto-optical rotation transit effect as a function of 
excited-state hyperfine splitting. Plotted are the component 
due to polarization produced by depopulation (dash-dotted) 
due to polarization produced by repopulation (dashed) and 
the total signal (solid). Parameters as in Fig. 14 
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i.e., linear in A e , with a modified dispersive shape that 
falls off far from resonance as 1/Af, where Ax is the 
detuning from the center of the F g = 1 — > F e transition 
group. 

The previous discussion also explains why the two 
peaks in the component due to depopulation seem to can- 
cel as they overlap, even though they have the same sign: 
the factors in the signal due to the creation and detec- 
tion of alignment cancel individually (Fig. |13| ; it is only 
in their product that the two peaks have the same sign. 

We now consider the case in which both the ground- 
and excited-state hyperfine splittings are small, so that 
all of the hfs is unresolved. To lowest order in A g and A e 
we have 



A, 



*A 



k sXi (T/ 2) 4 A 2 
24 [(r/2) 2 



(30) 



where A is the light detuning from the line center of the 
Dl transition. As we expect, the component of the signal 
due to polarization produced by repopulation is propor- 
tional to A g for small hyperfine splitting. The component 
of the signal resulting from depopulation-induced polar- 
ization also enters at this order. The optical rotation 
spectrum in this case is double-peaked, and falls off as 
1/A 6 (Fig. We I. 



B. Doppler-broadened transit effect 

We now consider an atomic ensemble with a 
Maxwellian velocity distribution, but a low rate of 
velocity-changing collisions, so that the atomic veloci- 
ties do not change between optical pumping and probing. 
This is the case for an atomic beam experiment, or for the 
"transit effect" in a dilute-vapor cell. Because the atoms 



FIG. 16: As Fig. |14| with A g and A e varied simultaneously 
(A g = 20A e ). 



have a fixed velocity, the signal from each velocity group 
can be found individually and then summed to find the 
total signal. Thus the signal from the Doppler-broadened 
transit effect is found by multiplying the Doppler-free 
signal found in the previous subsection by a Gaussian 
weighting function representing the Doppler distribution 
along the light propagation direction and then integrat- 
ing over atomic velocity. We can perform this integral 
analytically in different limiting cases. 

We first consider the commonly encountered exper- 
imental case in which the hyperfine splitting is much 
greater then the natural line width of the excited state, 
i.e., the Doppler-free spectrum is well resolved. In this 
case, for a given light frequency and atomic velocity, 
the light acts on at most one transition between hyper- 
fine levels. Thus the excited-state hyperfine coherences 
can be neglected, and the cancelation effects due to the 
overlap of resonance lines do not appear. As found in 
Eq. (A24l, the Doppler-free rotation spectrum then ap- 



pears as a collection of peaks, one centered at each op- 
tical resonance frequency, each with line-shape function 
f{oj'p cF ) = L(u> FsFs ) 2 , i.e., the square of a Lorentzian 
line shape. (One Lorentzian factor is due to optical 
pumping, the other to probing.) 

In this case, the Doppler-broadened signal is found by 
making the replacement / — > Jdb, where the velocity 
integral for fuB takes the form 

f DB (A FeFg ) = J dv k f(-A FcFg + k B v k )G(v k ), (31) 
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where 



G(v k ) 



c -(k B v k /r D ) 2 



(32) 



is the normalized distribution of atomic velocities along 
the light propagation direction k, ks is the Boltzmann 
constant, and Yd is the Dopplcr width. This integral can 
be evaluated in terms of the error function. Under the 
assumption r C To that we will employ here, the inte- 
gral can be approximated by replacing / with a properly 
normalized delta function, resulting in 



f DB (A FiiF ) Ri^^-e" 

4 1 D 



(A FeF /r D f 



(33) 



The Doppler-broadened spectrum, given explicitly by 
Eq. (A25l, thus consists of a collection of resonances, 
each with Gaussian line shape. For the Dl line with 
7 = 1/2, we have 



dip K S X 



576 



i f (3 + R) e -^W + 3e -(^/r D f\ 
i V / 



(34) 

Here £o is the absorption length for the Doppler- 
broadened case, given by 



4y^f T D (2 J g 
Rn\ 2 Y (2J e 



_1) 
1)' 



(35) 



Equation (|34|) is valid for A e ,A g ,T D > Y. Note that all 



the terms in this expression have the same sign; thus no 
cancelation occurs when the resonances overlap. This is 
because the Doppler-free resonances all have the same 
sign when the Doppler-free spectrum is well resolved 
(Fig. 14 1, so when the Doppler-broadened spectrum sam- 



ples more than one resonance, the contributions from 
each resonance add. 

The same approach can be generalized to describe the 
case in which some or all of the hyperfine splittings are on 
the order of or smaller than Y. In this case, the Doppler- 
free spectrum is not composed entirely of peaks with a 
shape given by f(u>' F F ). Nevertheless, as long as each 
resonance or group of resonances has frequency extent 
much less than the Doppler width, we can approximate 
it as a delta function times a coefficient given by the 
integral of the Doppler-free spectrum over the resonance. 
For the Dl line with I = 1/2 and ^4 e ,r < T D < A g , 
this procedure yields [Eq. ( |A26 1] 



dip _ Aln sXl (6 + i?)e~ A i. l/r ' 



576 (r 2 + A%) 



(36) 
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FIG. 17: Maximum of the spectrum of the Doppler-broadened 
effect as a function of excited-state hyperfine splitting. Plot- 
ted are the component due to polarization produced by de- 
population (dash-dotted), due to polarization produced by 
repopulation (dashed), and the total signal (solid). Equation 
p6| is used for A e < 10, and Eq. {34} is used for A e > 10. Pa- 
rameter values in units of P are To = 100, 7< 1, j4 s > Td- 



signal over the entire range of A e to excellent approxi 
mation, as verified by a numerical calculation. Figure 17 
shows the maximum of the rotation spectrum as a func- 
tion of the excited-state hyperfine splitting. As discussed 
above, as A e is reduced, there is no suppression of the 
optical rotation signal when the Doppler-broadened hfs 
becomes unresolved. Only when the Doppler-free spec- 
trum for a particular velocity group becomes unresolved 
is there suppression, as described in the previous subsec- 
tion. 

Spectra for the Doppler-broadened transit effect are 
shown in Fig 



18 



and for A g and /. 



for large A g and various values of A e 



varied together in Fig. 19 
In the case in which both A 
Doppler-free rotation spectrum is entirely of the same 
sign [see Eq. (30 1 and the bottom plot of Fig. 16 . The 



and A g are small, the 



Doppler-broadened signal thus behaves similarly to the 
Doppler-free signal, because no additional cancelation 
takes place upon integrating over the velocity distribu- 
tion. The signal for the Dl line with 7=1/2 and 
A g , A e , T < T D is given by 



dip 

] 1H 



A R A 



*A 



K s%i -A 2 /r 2 D 

96r 2e 



(37) 



Wall effect 



The rotation in this case goes as A\ for small A e ; the 
term linear in A e [Eq. (29 1] is odd in detuning and con- 
sequently cancels in the velocity integral. 



Since Eq. ( 34 ) applies when A e » T and Eq. ( 36 ) ap- 
plies when A e <fj), we have that — if Td is sufficiently 
larger than T — the two formulas together describe the 



We now consider systems in which the atomic veloc- 
ities change in between optical pumping and probing. 
This is the case for the "wall effect" in antirelaxation- 
coated vapor cells: atoms are optically pumped as they 
pass through the light beam, and then retain their polar- 
ization through many collisions with the cell walls before 
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FIG. 18: Spectra, as Fig. |14[ but for the Doppler-broadened 
transit effect. Parameter values in units of V are Td ~ 100, 
7< 1, A g > r D . 
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Since we are now describing the average over all of the 
atoms in the cell, and not just the illuminated region of 
the cell, we take 7 to be the average ground-state relax- 
ation rate for an atom in the cell, rather than the transit 
rate through the light beam. We also multiply the po- 
larization by the illuminated fraction of the cell volume, 
Miium./^coii (assuming this fraction is small), to account 
for the fact that the light pumps only some of the atoms 
at a time. 

For the specific case of the D 1 line for an atom with 
1=1/2, Eq. ( |A10| takes the form, given in Eq. of 
a linear combination of Lorentzian functions L(uPp e p ). 
This simple form arises because, due to the selection 
rules for this transition, no coherences are formed be- 
tween excited-state hyperfine levels. For a general sys- 
tem this is not the case; however, if the excited-state 
hyperfine splitting is greater than T, the excited-state 
hyperfine coherences are suppressed, and all resonances 
once again have Lorentzian line shapes. Thus, assuming 
that r <C Td, the velocity integral can be accomplished 
by replacing L(u' F F ) by 



r 

dv k L(-A FisF + kv k )G(v k ) w ——e 



-(A FeFg /T D f 



The polarization in this case is given by [Eq. (A27l] 



(38) 



R r 
3 



e (A li0 /r D ) 2 _ e -(A M /r D ) : 



(39) 



where the saturation parameter for the wall effect is de- 
fined by 



Hum. 



(40) 



We make this new definition because, in the wall effect, 
light of a single frequency illuminating just part of the 
cell effectively pumps all velocity groups in the entire cell. 

The signal due to each velocity group is given in terms 
of p (20) (l) by Eq. p6j; integrating over velocity to find 
the total signal, we obtain [Eq. (A28)] 



FIG. 19: As Fig. |18| with A g and A e varied simultaneously 
{A g = 20A e ). 



returning to the beam and being probed (Fig. 111. A 
similar situation occurs in vapor cells with buffer gas. 

We assume that the atomic velocities are completely 
randomized after optical pumping. Then the density ma- 
trix for each velocity group is the same; to lowest order 
in light power, we can find the velocity-averaged polar- 
ization by integrating the perturbative expression (A10) 



over velocity with the Gaussian weighting factor ( 32 1 



dtp 



3 r 



-(A ,i/r D ) 2 



-(A!,i/r D ) : 



*iP (20) (ll)- 



(41) 

The spectrum of the signal due to the wall effect is quite 
different than the spectrum of the Doppler-broadened 
transit effect signal, and is in a sense more similar to that 
of the Doppler-free transit effect [37] . Equations ( 39 1 and 
(Ell have the same form as the Doppler-free Eqs. (24 1 



and ( 26 ) , with Lorentzians of width T replaced by Gaus- 
sians of width Tr>. Thus the rotation signal produced 
by the wall effect has similar spectra and dependence on 
hyperfine splitting as the Doppler-free transit effect, but 
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FIG. 21: As Fig. |15| b ut for the wall effect. Parameters are 
the same as in Fig. |17| 



with scale set by the Doppler width rather than the nat- 
ural width. This is illustrated in Figs. [20] and |2"T] for the 
case of large ground-state hyperfine splitting. Figure [20] 
shows the optical rotation spectrum for various values of 
A e , and Fig. [21] shows the maximum of the rotation spec- 
trum as a function of A e . These figures can be compared 
to Figs. [14] and [T5] for the Doppler-free transit effect. In 
particular, we see the same phenomenon of two resonance 
peaks of the same sign appearing to cancel as they over- 
lap (observation of this effect in anti-relaxation-coated 
vapor cells is discussed in Ref. [37] and in buffer-gas cells 
in Ref. [35]). The explanation for this is the same as in 
the Doppler-free case. Also as in the Doppler-free case, 
the rotation is linear in A e to lowest order, and this lin- 
ear term is due to polarization produced by spontaneous 



dip 
] ~dl 



288 s lH r% 6 



(42) 



Spectra for the case in which A e and A g are varied to- 
gether are shown in Fig. [22] and are also similar to the 
Doppler-free transit effect (Fig. 16). When both A e and 
A g are small, the signal to lowest order in these quantities 
is given by 



R 

3 9 



K S X\A 2 



-2A 2 /rj, 



(43) 



D. Higher nuclear spin and the D2 line 

When nuclear spins I > 1/2 are considered, several 
complications arise. The clearest of these is that the two 
ground states now have angular momenta F g — I±\ > 1, 
so that they can both support atomic alignment and pro- 
duce optical rotation. A more subtle difference is that, 
with higher angular momenta in the excited state, co- 
herences between the excited state hyperfine levels can 
be created when the excited-state hyperfine splitting is 
on the order of the natural width or smaller. (Ground- 
state hyperfine coherences can be neglected as long as the 
ground-state hyperfine splitting is much larger than the 
ground-state relaxation rate.) This can change the opti- 
cal rotation spectrum, and also causes the symmetry be- 
tween the Doppler-free transit and wall effects discussed 
above to be partially broken, as we see below. 

However, many of the results obtained above for the 
1=1/2 system are a consequence of the general argu- 
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ments discussed in Sec. [TT] and thus hold for any nuclear 
spin. In particular, the dependence of the optical rota- 
tion signal on the hyperfine splitting for large ground- 
state and small excited state splitting [Eqs. (29), (36 1, 
and ( 42 1] and for both ground- and excited-state hyper- 
fine splitting small [Eqs. (30 1, (37 1, and (43 1] remains the 
same. We have, for large Ag and small A e , and for a par- 
ticular transition group, the following three expressions. 
For the Doppler-free transit effect: 



, (r/2) 4 A Fs 
to— tx A e K s x F R- 



[(r/2) 2 

for the Doppler-broadened transit effect: 



A 2 



(44) 



, -a% /r 2 D 



and for the wall effect: 



dL P A T} Ap o e 

~dt A e K sXp s R ! ^2 



2(A F jr D f 



(45) 



(46) 



For A g and A e both small, we have, for the Doppler-free 
transit effect: 



t ^ oc A P | .1, 



-A g 
3 9 



« s (r/2) 4 A 2 



[(r/2) 2 

for the Doppler-broadened effect: 



(47) 



R A \ K e 

3 V f^ e 



and for the wall effect: 



o dL P A 
t — cx A e 



A R 



3 9 



t s A 2 

r 4 

1 D 



,-2A 2 /r 2 D 



(49) 



To illustrate the differences that arise when the nu- 



clear spin is increased, we plot (analogously to Figs. 15 
17 and 21 1 in Fig. 23 the maximum of the rotation spec- 
tra for large A g as a function of A e for the Doppler- 
free transit, Doppler-broadened transit, and wall effects. 
Three values of the nuclear spin are used, I = 1/2, 3/2, 
and 5/2, and for 1 = 3/2 and 5/2 the rotation on the 
F g — I ± 1/2 — + F e lines is plotted separately. Rotation 
due to polarization produced by the depopulation and 
repopulation mechanisms is plotted, as well as the total 
rotation signal. In many cases these two contributions 
are of opposite sign, so the details of the total signal can 
depend on how closely the two contributions cancel each 
other. (The cancelation tends to be more complete for 
the F g = I — 1/2 lines.) However, the qualitative features 
of these plots follow, in large part, the pattern exhibited 
in the 7=1/2 case. One exception is the behavior of the 
wall effect plot for A e in the neighborhood of the natural 
width. As mentioned above, when I > 1/2, excited-state 



hyperfine coherences can form when the excited-state hy- 
perfine splitting becomes small. This leads to "interfer- 
ence" effects when the Doppler-free resonance lines over- 
lap that do not occur when the Doppler-broadened res- 
onance lines in the wall effect overlap. This breaks the 
symmetry between the wall effect and the Doppler-free 
transit effect that is found in the 7=1/2 case. 

We now discuss the J g = 1/2 —> J e = 3/2 D2 transi- 
tion. The presence of three hyperfine levels in the excited 
state leads to additional features in the dependence of the 
signal on the hyperfine splitting (Fig. 24 1. However, the 
fact that the ground-state electronic momentum is still 
J g = 1/2 means that the dependence of the signal on 
the excited-state hyperfine splitting as A e goes to zero 
remains the same, for the reasons discussed in Sec. [IT] 
Thus, to lowest order in A e , the rotation signals on the 
D2 line for large A g are given by Eqs. (44 46 1. (We set 



the hyperfine coefficient B e to zero for simplicity.) 

Considering the signals obtained when both the 
excited- and ground-state hyperfine splittings are small, 
we expect somewhat different behavior for the contribu- 
tion due to polarization produced by repopulation pump- 
ing than in the Dl case. This is because the excited- 
state electronic angular momentum is J e = 3/2, so that 
production of rank k = 2 < 2J e atomic alignment in 
the ground-state by spontaneous emission is allowed even 
when the ground-state hfs is unresolved (Sec. II B I. The 
lowest order dependence on hyperfine splitting for the D2 
line is given by 



4 -jr oc A e 
at 



A R -R 2A 



3 9 



(50) 



for each of the three effects, with the spectral line shapes 
remaining as in Eqs. (47j49|. Note that there is now a 
term that depends on polarization due to repopulation 
that does not go to zero as A g goes to zero. 



E. The alkalis 

We now examine the consequences of the preceding dis- 
cussion for the alkali atoms commonly used in nonlinear 
magneto-optical experiments. In Fig.[25]the maximum of 
the spectrum of optical rotation is plotted for the Dl and 
D2 lines of several alkali atoms. The Doppler-broadened 
transit effect is shown in Fig. |25~H and the wall effect is 
shown in Fig. 



25 }. (Numerical convolution was used to 



obtain these results, because the alkalis do not all satisfy 
the conditions under which the analytic formulas were 
derived.) The nuclear spins, hyperfine splittings, excited 
state lifetimes, and Doppler widths all vary between the 
different alkali atoms. However, focusing our attention 
on the hyperfine splittings, which have the greatest de- 
gree of variation, we can see the correspondence of these 
results to the preceding discussion. In particular, we have 
seen that the magnitude of the Doppler-broadened tran- 
sit effect is largely independent of the hyperfine splitting 
when the splittings are greater than the natural width of 
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FIG. 23: Maximum of the normalized optical rotation spectra £o/(K s XF g )dtp/dl for the Doppler-free transit, Doppler-broadened 
transit, and wall effects on the Dl line for I = 1/2, 3/2, and 5/2. We assume Td ~ 100, 7 C 1, A s > Fd in units of V. The 
maxima for the F g = I ± 1/2 — » F e transitions are plotted separately. Each plot shows rotation due to polarization produced 
by depopulation (dot-dashed line), rotation due to polarization produced by repopulation (dashed line), and total rotation 
(solid line). The sharp features seen in some of the plots occur when two terms contributing to the largest resonance cancel. 
In general, since all resonances are not canceled at the same time, the maximum of the spectrum does not go to zero. 



the transition. This is generally the case for the alkalis, 
leading to the relative constancy of the magnitude of the 
transit effect among the alkalis. For the wall effect, on 
the other hand, we have found that the magnitude of the 
effect diminishes when the hyperfine splitting becomes 
less than the Doppler width. In the alkalis the excited 
state hyperfine splitting is generally on the order of or 
smaller than the Doppler width, and the general trend 
is that the ratio of hyperfine splitting to Doppler width 
increases as the atomic mass number increases. This ac- 



counts for the general upward trend in Fig. |25b . The 
trend is not completely consistent: the hyperfine split- 
ting of K is smaller than that of Na, which is reflected in 
the plot of the wall effect. 



IV. CONCLUSION 

In experiments involving light-induced polarization in 
the alkali atoms, the effect of partially resolved hyper- 
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FIG. 24: Maximum of the normalized optical rotation spectra, as in Fig. 23 but for the D2 line. For the I = 3/2, F g = 2 and 
the I = 5/2, F g = 3 systems for the Doppler-broadened transit effect, the two contributions to optical rotation nearly cancel, 
with the consequence that the approximations used in obtaining the analytic formulas for the total Doppler-broadened signal 
begin to break down. Numerical convolution is employed in these cases. 



fine structure is of practical importance. We have ad- 
dressed this question from both descriptive and quanti- 
tative standpoints. We have formulated rules describing 
various restrictions on the rank of atomic polarization 
moments that can be created or detected by light in cases 
when cither the ground- or excited-state hyperfine struc- 
ture is completely unresolved. We have also studied the 
particular situation of nonlinear Faraday rotation under 
various experimental conditions in more generality, and 
presented analytic formulas giving the results of optical 
rotation measurements when the hfs is unresolved, par- 
tially resolved, or completely resolved. 
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23 Na 39 K 85 Rb 87 Rb 133 Cs 23 Na 39 K 85 Rb 87 Rb 133 Cs 

FIG. 25: Maximum of the spectrum of normalized optical rotation for the (a) Doppler-broadened transit effect and (b) wall 
effect for various alkali atoms. Circles indicate the Dl line and triangles indicate the D2 line. Room temperature Maxwellian 
velocity distributions are assumed. Normalized rotation is defined here as £o/(K s xi+j g )(dtp/d£), where lo in this case is the 
unsaturated absorption length at the detuning that gives maximum absorption. The normalized magnitude of unsuppressed 
optical rotation is nominally on the order of unity; however, this is to some degree dependent on the normalization convention 
chosen. For example, if the maximum matrix element of d z is used in the definition of k 3 , rather than the reduced matrix 
element, the values in this plot are increased by a factor of ~6. 



APPENDIX A: NONLINEAR 
MAGNETO-OPTICAL ROTATION WITH 
HYPERFINE STRUCTURE 

1. Perturbation theory with polarization moments 

The time evolution of the atomic density matrix p un- 
der the action of a time-independent Hamiltonian H is 
given by the Liouville equation, which can be derived 
from the Schrodinger equation (with phenomenological 
relaxation terms added by hand). Setting h = c = 1, the 
Liouville equation can be written 



1 



P : 



i[H,p]--{T,p} + A + Tr(Fp) 



(Al) 



where the square brackets denote the commutator and 
the curly brackets the anticommutator, L is a Hermitian 
relaxation matrix that accounts for relaxation mecha- 
nisms such as transit relaxation due to atoms leaving the 
region of interest and intrinsic relaxation of excited states 
due to spontaneous emission, A accounts for repopulation 
mechanisms, such as transit repopulation, that do not de- 
pend on p, and T is the spontaneous emission operator, 
accounting for repopulation of ground states due to spon- 
taneous emission from excited states. (We neglect other 
relaxation and repopulation mechanisms, such as spin- 
exchanging collisions, which may require the inclusion of 
additional terms). The spontaneous emission operator, 
denned by [30JET] 



(A2) 



connects a pair of excited states |s), \r) to a pair of ground 
states |m), |n); the trace in Eq. Al is taken over the 
excited state pair. 

We can expand the operators appearing in the Li- 
ouville equation in terms of the polarization operators 
T^\f 1 F 2 ) according to 



4 = ^AW(F 1 F 2 )TW(ir 1 F 2 ), 



(A3) 



where % runs over all pairs of states in the system. 
(Here F is understood to represent the total angular mo- 
mentum quantum number as well as any additional quan- 
tum numbers necessary to distinguish between two states 
with the same total angular momentum.) The expansion 
coefficients A^^Fq) can be found from the Wigner- 
Eckart theorem, along with the orthonormality condition 



Tr(TW(^F 2 ) (t(/\f{F£ 
and the phase convention 



3kk' Sqq' Sp 1 F[ $F 2 F£ 

(A4) 



T^{F 1 F 2 )) = (-1) F ^ F ^T^(F 2 F 1 ) 
= T {Kq \F x F 2 ). 



(A5) 



There are several equivalent expressions for the expansion 
coefficients; one such form is 
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A(^)(F 1 F 2 ) = ^(-1) 



F 1 -F 2 +q 



2k + 1 
2F 2 + 1 



(FxmK, -q\F 2 m')A Fim . F2 . 



(A6) 



The set of expansion coefficients for the density matrix are known as polarization moments. Performing the 
expansion of each operator, and using appropriate tensor product and sum rules, the equation of motion for the 
polarization moments is found from the Liouville equation to be 



p^ q \FtF 2 ) = - ^ y/(2 K ' + 1)(2k" + l)(K'q' K "q"\Kq) 



ff("V)( FlF 3) _ l -T^^(F 1 F 3 ))p^"""HF 3 F 2 ) 



(FtFa) [H^"""\f 3 F 2 ) + l -T^"""\F 3 F 2 ) 



(A7) 



+ A™{F 1 F a ) + | W §X;(^l|d|l^)^ ) M)(^||d||F 2 )(-l)^ + '»+ 1 |* 2 f'\ ' 

where all variables not appearing on the left-hand side are summed over (the variables F e and F' e appearing in the 
last term relate to spontaneous emission and run over only those states of higher energy than F\ 2 )- Here the arrays 
enclosed in in curly brackets are the 6j symbols. 

We now suppose that the total Hamiltonian H = H + V, where H is diagonal and V is a time-independent 
perturbation. We also assume that T and A are diagonal. More precisely, we assume that only T^ 00 ^(FF), A( 00 \FF), 
and H(°°\FF) are nonzero (for arbitrary F). Taking the steady-state limit in Eq. ( |A7| and expanding to second 
order in the perturbation V, we find for a ground-state polarization moment 



p^\F g F g ) = 



7 



lLV FF N q 



6 K0 5 q0 V^Fi+T - (_i)2f s +«'+-" A /(2 K ' + l)(2K"+l)( K 'g , KVI«g) 

x W K ' K \ u Fg F>+u F , Fg (wV) ((tV0 ^ 

[F g F g F'j u FgF fi FgF ,u F , Fg 

4 cj (-1) 2F i' +2F 2+«+«'+«" v /(2k' + 1)(2k" 4 
3 

x (F 2 ||rf||F g )(F g ||d||F() ~ F (re " 9 " ) (F 1 F / )F( K ' 3 ' ) (F / ^) 



f/e" 


«/ 


.1 


w 


'k F„ 


*.\ 


1*5 


Fi 






[1 "{ 





(A8) 



Here we have neglected the possibility of cascade decays 
and assumed that V does not couple a state to itself. We 
also have assumed that A repopulates all ground-state 
sublevels equally; i.e. A Fm>Fm — j/N g , where 7 is the 
ground-state relaxation rate and N g is the total number 
of ground-state sublevels. The complex energy splitting 
uj FiF2 is given by 

5 FlF2 = E Fx -E F2 -~ (T Fl + T F2 ) , (A9) 

where E F = (2F + 1)-^ 2 H^ 00) (FF ) is the unperturbed 
energy and T F = (2F + l)-V2r(°o)(FF) is the total 
relaxation rate of a state F. 



2. Doppler-free transit effect 



We now apply the results obtained in Sec 
the three-stage calculation described in Sec 



Al 



III A 



to 
In 

stage (a), we consider a z-polarized optical electric field 
S = SoRe(ee l ' k ' r ~ UJ *'), where £q is the optical electric 
field amplitude, e = z is the polarization vector, k = x is 
the wave vector, and uo is the optical frequency. We let V 
represent the electric-dipole Hamiltonian in the rotating- 
wave approximation: V = — ^d z £o- (Here the prime 
refers to the rotating frame.) We assume that the mag- 
netic field is absent in this stage. From Eq. (A8) we find 
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£(-*) 



f,-f. k (2F e + l)(2F g + l) / I+2Ja 
Ks (2/ + l)(2J. + l) ( ij 



R E (-l) F ^^(2Je + 1)(2F; + 1)(2^ + 1) j 



1 1 2 

v F g F e 

112] 

F P F> F'\ 



J e F e I 
F g J g 1 



F 9 F g 2 
F P FL 1 



L (u F v) 



J, F, i] \j e F e i] fx FL A \j e FL i] L{u' FaF ,)L{u' F , F X 



F g J g 1 



F' J 1 

r g J g 1 



F g J g 1 



Fg J g 1J i( F'^F'F' 



(A10) 



where all variables are as defined in Sec. IIII Al We 
have evaluated matrix elements using the Wigner-Eckart 
theorem and have used the relation (see, for example, 
Ref. [35]) 



where K IJF = F(F + 1) - 1(1 + 1) - J(J + 1) and Aj 
and Bj are the hyperfine coefficients. The last term is 
zero for J < 1/2 or I < 1/2. 



RT = 



{Jg\\d\\J e f 



3 2 J e + 1 

The unperturbed energies can be evaluated with 

Ejfm = Ej + -Kij F Aj 

3 K IJF (K IJF + 1) - § J(J + 1) J( J + 1) 



(All) 



1(21 -l)J(2J - 1) 



5 j, 
{A43} 



In the case in which the excited-state hfs is well 
resolved in the Doppler-free spectrum (u>f c F' 3> T), 
Eq. (|ATo|) reduces to 



+ Y^R (-l)^- Fc (2J e + l)(2i^ + 1)(2F C + 1)1 

F' I 



1 1 2 II F g F g 2 II J e F e 1 1 I J e F e I 

Fe F e F'g J \F e F e 1 \F g J g 1 \F> J g 1 



L(^ FcF ,)j. 
(A13) 



In stage (b), the ground-state density matrix, which is 

initially in the state found in stage (a) , evolves under the 

influence of a magnetic field _Bx. We will require only 

(21) 

the value of the polarization moment p b _(F g F g ). Using 
the Hamiltonian Hg = — fi ■ B in Eq. ( |A7[ ) and solving 
for the steady state, we find 



P ( b 21) (F g F g ) 



''2V2 



x F Jr\F g F g ). 



(A14) 



the wave equation: 



d£_ 



d£_ 



47T 



dP 



(A15) 



where £ is the distance along the light propagation direc- 
tion, and P = n Tr pd is the medium dipole polarization 
(n is atomic density), which can be written in terms of 
the rotating-frame density matrix p' as 



where the magnetic-resonance line-shape parameter XF g 
is defined in Eq. (27 1. 



In stage (c) the ground-state polarization is probed. 
The effect of the atoms on the light polarization as the 
light traverses the atomic medium can be found in terms 
of coherences between ground and excited states using 



P = n 2 Re 



H „i(k-r-wt) 



pm^-rnp 



(A16) 



nip 



where m runs over ground states, and p runs over excited 
states. Using the parameterization of a general optical 
electric field in terms of the polarization angle if and 
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ellipticity e, 

S = £ Re (e^-r-^+<f>) [( cos 99 cos e - i sin ip sin e) e x 
+ (sin </? cos e + i cos ^ sin e) §2] 



(A17) 

where ei^ are orthogonal transverse unit vectors, we ob- 
tain for optical rotation in the case of linear polarization 



dip 

lie 



So 



E Im Ppm d mp • (k x e) 



(A18) 



Using first-order perturbation theory for the optical co- 
herences and neglecting coherences between nondegener- 



ate ground states, we obtain the optical rotation for weak 
probe light in terms of the ground-state density matrix: 



dip 



-27ro;nIm 
where we have defined 



f j =J2 



e • (3 • (k x e) 



i ■■ : ! > : ■: jP m 



(A19) 



(A20) 



Expanding the tensor /3 in terms of the ground-state po- 
larization moments, we obtain 



&= E e- q ,i- q „(lq'lq"\K,</ + q") 



F g F eK q'q" 



J F r F a 



1 1 K 

Fg F e 



\(F g \\d\\Fe)\ 2 P (K ' q ' +q "HF g F g ), (A21) 



where e„ are the spherical basis vectors. 



Evaluating (A19) for our case and using Eq. (A14) gives 



- E M) Fs+F w + i)m + i)(2j 9 + 1) I 



1 1 2 I I J e F e I 

Fg Fg fJ] F n J„ 1 



L(u') XF j™\FgFg), 



g "9 



(A22) 



where the unsaturated absorption length for the J g — > J e transition is defined in Eq. (28 1. Substituting in Eq. (AlOl 
results in the complete expression for optical rotation for the Doppler-free transit effect: 



to 



dip 3 
~dl = 4 



*s E 

F g F e F{J 



2F 0+ f"^ (2^e + mF'J + l)(2F g + l) 2 f 1 1 2 If J e F e " J 



(2/ + 1) 



1 1 2 I I J e F e I 

Fg Fg F e I I F„ J fl 1 



Fq I I -f O ^(7 1 



XF a 



L {u FeF )L(Jp,, F ) 



g °g 



R £ (-l)^-^ (2 J e + 1)(2F; + 1)(2^ + 1) j 



1 1 2 \\Fg Fg 2 

F e F' e F'J\F e F' e 1 



j e F e A \j e f„ 1} fx fl 1] \ j e f: 1} Hu Fc F>) L (u F >F>) L (u' F ,, F y 



Fq J a 1 I I Fg Jg 1 I I Fq J q 1 I I F n Jg 1 



g °g 



(A23) 



For completely resolved hfs (ujp^pi ,uip pi ^> T), this reduces to 



(2F e + l) 3 (2F 3 + l) 3 I 1 1 2 II J e F e I 



F a F e 



(2J+1) 



XF a 



g u g 



(-1) 



2/+2J„ 



I 1 12 

(2F e + l)(2F 5 + l) F g F e 



+ R(-l)*'-*-(2J e + l) 



1 1 2 \\F g F g 2\\j e F e I 

Fe F e Fg \F e F e 1 \Fg Jg 1! 



(A24) 
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3. Doppler-broadened transit effect 



The procedure used to obtain the optical rotation signal in the Doppler-broadened case is described in 
Sec. |IIIB| When the ground- and excited-state hypcrfinc splittings are al l mu ch greater than the natural width 
(u>P e p^,u!p g p/ ,Tu 3> r) we have, applying the integration procedure to Eq. (A24l, 



3 v Fo (2F e + l) 3 (2F 3 + l) 3 f 1 1 2 lf 



Je F e I 

F 9 J 9 1 



e -(A FcFg /r D y XF 



(-1) 



2/+2,/„ 



1 1 



+ i?(-l) F ^(2J e + l) 



(2F e + l)(2F g + l) }F q F q F e 



1 1 2 

F& Fp_ Fa 



F 9 F g 2 
F„ F K 1 



Je F e I I 
^ 4, 1J 



(A25) 



where the unsaturated absorption length for the Doppler-broadened case is given by Eq. ( 35 ) . 

In another limit in which the ground-state hyperfine splittings are much greater than the natural width and the 
excited-state splittings are much less than the Doppler width (u>F e p^,T <C Tr>, T <C ^F g F')i we have 



dip 3 
] ~d£ = 8 



«. E (-d 2F! 



+F ,j^ Fc (2F e + l)(2Fj' + l)(2F g + l) 2 J 1 1 2 J e f e / J e / 



F„F,F'J 



x (-1) 



F F F" 
r g r g r e 



2/+2J„ 



1 1 2 

F g F g F ( 



}+i?E(-l) 



(27+1) 

' - F =(2J e + l)(2F g + l)(2^' + l) 



F g J g 1 



F g J g 1 



X F C , 



F g J g 1 



^2r 4 + (2r 2 + U FeF ,)u Fe F»WF£FZ 



2(r 2 + c4 F ,)(r 2 + WF , F „) 



1 1 2 

F fi F' F n 



-(A Fg /r D ) 2 r 2 

r 2 + c4 F „ • 



^ 9 F g 2 
F R F' 1 



(A26) 



4. Wall effect 



The procedure for obtaining the signal in the wall effect case is described in Sec. IIIC For excited-state hyperfine 
splittings much greater than the natural width (T <C u>f c f' ,^d), we have for the ground-state polarization 



+ i?^(-l) F ^(2J e + l)(2F g + l)(2F e + 1) J 



-(A FeFg /r D y 



1 1 2 F„ F„ 2 \ J e F K I 



a ^ a 

F e F e F' g \ [F e F e 11 IF; J g 1 



(A27) 



where the saturation parameter for the wall effect is defined by Eq. (40 1. The optical rotation signal is then given by 

\ f \2 



^ = -^E(-l)^(2F fl + l)(2Fi + l)(2J fl + l)|^ ^ |,J|^ J J| e-^/r°r XFgP m {Fg F g ) 



dt 



3v^r 



V- ( U 2F + F'-F C ( 2F e + \){2F' e + l)(2F g + l) 2 / 1 1 2l jj e F e iYfje K l\ 

( 2/+1 ) F 2 U Ja 1/ U ^ 1/ 



x z Fs e- (A W rD)2 (-1) 2 '+ 2 ^ 112 e-( A -e^/r D ) 2 



1 1 2 M F s F g 2 M J e F e J I e -(A FcF ,/r D ) 2 
F e F e F' f 1 F R F„ 1 [ ] F' J a 1 I 



(A28) 
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